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ABSTRACT 

Ion acoustic waves have been artificially excited and detected in 
the Steady State Plasma Facility at the Naval Postgraduate School by a 
magnetic excitation device and an optical detection system. The dis- 
persion of these waves has been studied as a function of the exciting 
frequency and the confining static magnetic field. 

Dispersion curves are presented which compare the experimental data 
to a new dispersion equation which is derived in the present thesis. 
After the validity of this dispersion relation is demonstrated, the 
measurements are used as a diagnostic tool, and five operating para- 
meters of the plasma facility are deduced from the measurements. 

Theories are presented to explain the behavior of the dispersion 


curves and three of the five operating parameters. 
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Section 1 
INTRODUCTION 

For the purposes of the present thesis a plasma may be defined as 
a quasi-neutral electrically conducting mixture of positive ions and 
electrons plus a number of neutral background gas atoms. Due to the 
particular plasma production mechanism employed in the Steady State 
Plasma Facility of the Naval Postgraduate School, the electrons are at 
a relatively high temperature (22 electron volts) on the cylindrical 
beam axis, while the ions have thermal energies of the order of .09 ev 
and the neutral gas atom energies are approximately .01 ev. 

Due to the high electrostatic restoring forces created by even a 
small internal charge imbalance, a plasma has the ability to maintain 
overall internal charge neutrality. If a charge imbalance should 
develop due to the random thermal motion of the electrons and ions, the 
particles will rearrange themselves to shield the remainder of the 
plasma from the resulting electric field. Thus beyond a certain average 
distance from the imbalance the electric field so created will not be 
felt. This distance is called the Debye length and is, of course, a 
function of the plasma parameters. In a true plasma the Debye length 
must be much less than the smallest dimension of the plasma system. 


The Debye length is given by: H 


where 
€. 15 the permittivity of free space 
K is the Boltzmann constant 


T is the electron temperature 


N is the equilibrium plasma density 
e 1s the electron charge 
Ion acoustic waves are longitudinal plasma density oscillations 

occurring at frequencies well below those of the electron plasma oscil- 
lations. The plasma ions and electrons move exactly in phase only at 
(9 - 0; for 42 # 0 in-phase motion is only approximately realized. For 
the frequencies employed in the present investigation, however, the ap- 
proximation is an excellent one. The ion oscillations are slightly 
larger in amplitude than those of the electrons. The ratio of ion per- 
turbation velocity to electron perturbation velocity can be shown to be 


1 + (Хр 
2ግኘ 





)2 B where ይ is the wavenumber. In Section 2 it will be 


2% 


Due to their small mass and high mobility the electrons see these 


shown that ( , so it may be neglected compared to 1. 





low frequency waves as quasi-static perturbations. Ton acoustic waves 
are Similar to ordinary sound waves in a neutral gas, except that the 
propagation mechanism is the Coulomb interaction between charged particles 
rather than the short range collision forces between molecules. Their 
existence was first predicted by Tonks and Langmuir in 1929 based on the 
fluid analysis of a two component plasma. Bg 
Before proceeding to an account of previous work, the Tonks- 

Langmuir dispersion relation for ion acoustic waves will be derived. 
The following assumptions are made: 

1) The plasma electrons have negligible momentum compared to 
the plasma ions. 

2) The plasma ions have zero temperature and are singly 
charged. Thus there is no ion pressure. 


3) There is always 100% ionization within the plasma. 
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4) There is no electric field associated with the wave. 
5) The Debye length is much smaller than the wavelength, 
that is እ ህ«-ሟ፦ : 
6) The drift velocity of the ions is zero. 
7) The unperturbed plasma has uniform density. 
Under these assumptions the momentum transfer equation for the ions can 


be written as: 





4 2 B - 
ZYP + р) = > Ша σ]]α d oy - (8 + 09 


where 

M is the ion mass 

p is the electron pressure perturbation 

P is the equilibrium electron pressure 

V is the ion velocity inducted by the wave 

N is the equilibrium plasma density of charged particles 
n is the plasma density perturbation of charged particles 
e is the electron charge 
56 is the static magnetic field strength 

z is a unit vector parallel to the direction of Be. 

Subtracting the momentum transfer equation for the ions in the un- 


perturbed system from the above equation and linearizing to the limit of 


small perturbations gives: 





nb q _ > - ^ ^ 
SNP = yr OV) -ANVgr *-Q.NV.9 (1) 
° е eBg 

where [l-is the ion cyclotron frequency = M 


Now for perturbations of the form ei bz σα. equation (1) can 


be solved to give: 
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в = М > r SE зру. ) 


апа 
NVg - app. “а? € ) 
Therefore 
N: (NVer - NV, 8)- < ы: рау — 








Taking the divergence of equation (1) gives: 


3% - -ieNL. (NV) - (QNS (Wer - NVQO), 


Or 


21 ב‎ = а? 22} 1 2p 1 
መ P = -іс)Ҹ. (МУ) - м S72 NE = pe 





(2) 


But by the linearized equation of continuity, 


NIIS) SE 3з- = 16дп 


Therefore equation (2) becomes: 


2. = wed 2, 4 -02 422, αἱ ο “1 | 
E vu бә "7 Sr ο a, 





Assuming there is no radial electron temperature variation, 








2 2 
J 2р = er qan, ὃ Б = זא‎ 2 n , >p = kr “Оп 
r 








dr? Ər Sr ' 
and 
а + : ብኪ. - Д?а) з -еУ п 5 = mn 
where 
с2 з E (definition) 


and T is the electron temperature. 
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Rearranging and defining 


2 Q2 /—— —À 


this equation becomes: 


Its solution is: 
me a, ci г Ф Маа 


where Jg is the zero order Bessel function. 


Therefore it must be true that 


= (2.405...)2 





2 
2 x. ρω... 
үк тт 
where the boundary condition that the perturbations vanish at the 


position r = / is used. (The use of this type of boundary condition 


will be more fully discussed in Section 2.) This can be rewritten as: 





2 5.79 Я -- 
: 
pd 426. 1) 
ma. 2. «У tt) 
Now for an infinite medium Де апа д; = = Or he (3) 


Thus for an infinite medium the application of a static magnetic field 
has no effect. Equation (3) is the well known Tonks-Langmuir dispersion 
relation for ion acoustic waves. 

From an elementary acoustic analogy the "hot" electrons provide the 
perturbing and restoring pressure gradients, while the massive ions pro- 


vide the inertia. This dispersion relation predicts a phase velocity of 
1 

ሬኃ = с = (KT). Also, since ο - дд ; the group velocity equals 

= M б oL 

the phase velocity and, there is no dispersion. 


A similar dispersion relation will be derived in greater detail 
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including additional effects in Section 2. 

Ion acoustic waves can be self-excited or externally excited both 
in gaseous discharges of low percentage ionization and in highly ionized 
plasmas. The existence of self-excited waves in low density gaseous 
discharges has been observed by Alexeff and Neidigh ῃ and by Little and 
Jones [5] in the form of standing waves in the discharge tubes employed. 
The present investigation deals exclusively with propagating waves. Two 
earlier experiments (one on standing waves and one on propagating waves) 
will be described. 

Ion acoustic waves have been excited by electrostatic methods (grids 
[6] and probes 7} and by magnetic perturbations of the local electron 
diffusion coefficient [5]. The magnetic perturbation method was used 
in this experiment. Its theory of operation appears in Section 3. 

Once the ion acoustic waves have been excited, they can be detected 
by grids, probes, and light sensitive detectors. In this experiment both 
probe detection and photomultiplier detection were used, although only 
the photomultiplier detection system was used to make quantitative mea- 
Surements. Detector theory and operation are also covered in Section 3. 

Wavelength and attenuation measurements were made as functions of 
exciting frequency for various static magnetic field strengths. The 
validity of the dispersion relation is demonstrated, the data are anal- 
yzed and the results presented as dispersion and attenuation curves. 

The new dispersion and attenuation theories derived in Section 2 
contain 5 adjustable parameters. Comparison of the experimental data 
with theory allows computation of these 5 parameters, whose values are 
presented as functions of the static magnetic field strength. The ion 


acoustic waves are thus employed as a diagnostic tool. Physical inter- 
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pretations of the experimental results (dispersion curves, attenuation 
curves, and parameter behavior) are also given in Section 5. 

Wong, Motley, and D'Angelo have carried out experiments on the dis- 
persion of ion acoustic waves in cesium and potassium plasmas in strong 
magnetic fields. [| 

Fried and Gould [з] have shown that the Vlasov equation predicts 
ion acoustic wave propagation even in the absence of collisions in the 
plasma. Under this condition, however, the analysis predicts that the 
waves will be strongly damped due to a "trapping" of ions moving with 
velocities close to the wave phase velocity. This can occur if there are 
ions whose temperatures are close to the electron temperature (T =y). 
From equation (3) the wave phase velocity is c = ску уЗ, and if T = Y 
there will be ions meeting this condition. 

This type of collisionless damping was first predicted by Landau. [9] 
It was observed by Wong, Motley, and D'Angelo, and their work provided 
the first experimental proof of its existence. The conditions of their 
experiment were such that ion-neütral collisions were negligible (the 
percentage ionization varied between 40% and 90%). Their plasma was 
produced by surface ionization of cesium and potassium atoms on a hot 
tungsten plate. This method guaranteed that electron and ion tempera- 
tures were approximately equal. The damping was so strong that standing 
waves were not observed. 

The Q-3 Machine (Princeton) was used for the experiment. The fea- 
tures essential to the experiment appear in Fig. 1. The plasma was pro- 
duced at one end of the machine, allowed to drift longitudinally, and 
terminated approximately 90 cm away by a similar tungsten plate. Radial 


confinement was achieved by adjusting the longitudinal magnetic field to 


15 


1111 
1 


1 ШІ ШІ | 







еә 380Ч4 9NIAI393U 


NV 38 
IA EES I 4 uo mA, ገሃህ1በ38 
М3169МП1 | | | 4רצפא+‎ | 


LOH == - - - - + - --—4--4------- 


es» 39084 NOILVLIIX3 ай 31V1d NA3LSONNL LOH 


πμ NN ; 2 мземунэ 
КЦ ШШ 









pts bm ሃበበጋሃሄለ 
[111 | mi 


TTE 

ІШ 
፥1 

$ 1111 | ШІПШІШ 








38084 רטאפאוחוט‎ QION 310S 





olv TUJSO 


“ЕГЕП! ЗМІНЗҮМ 6-0 30 ሸሃህ9ሃ!|!0 29ILVW3HOS 





І 913 


16 


12,000 gauss for cesium and 6000 gauss for potassium (these values were 
chosen to avoid possible ion cyclotron resonance). 

In the Princeton experiment an electrostatic excitation method was 
used, in constrast to the magnetic excitation method employed in the pre- 
sent investigation. The perturbing signal was applied to the excitation 
grid and modulated the percentage ion transmission through the grid. The 
movable receiving grid detected the density perturbations. 

Phase velocity calculations were made by measuring the slope of the 
curve of phase delay vs. receiving grid position (the movable grid) for 
the propagating wave. One result of the experiment was that the phase 
velocity was found to be independent of frequency for the range from 15 
ke to 100 ke (ion cyclotron frequency was 140 ke for cesium and 240 kc 
for potassium). This conforms to the predictions of the dispersion re- 
lation derived in the present thesis, although the frequency range rela- 
tive to the ion cyclotron frequency is different. 

The damping observed by Wong, Motley, and D'Angelo was found to be 
in excellent agreement with the theory of Landau damping. In Section 2 
the criteria necessary for the presence of Landau damping will be dis- 
cussed, and it will be shown that it can be ignored in the present inves- 
tigation. 

Figs. 2 and 3 show phase velocities for the ion acoustic waves in 
the Q-3 experiment for propagation parallel and anti-parallel to the plas- 
ma drift velocity. 

The drift velocity of the plasma was estimated from the difference 
between the phase velocities of the waves for parallel and anti-parallel 
propagation. In the present investigation the plasma drift velocity is 


one of the 5 adjustable parameters. The magnitude of the drift velocity 
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is machine calculated from the dispersion data. Both parailel and anti- 
parallel propagation at 1000 gauss were studied in the present investi- 
gation, and the difference in phase velocities for the two directions of 
propagation provided an independent determination of the drift velocity. 
The method is the same as that employed by Wong, Motley, and D'Angelo. 
This determination is discussed in Section 5. 

Little and Jones have made dispersion measurements on externally 
excited ion acoustic waves in the positive column of a iow pressure mer- 
сиг arc. B The percentage ionization in the discharge was approximately 
1%, the current through the mercury vapor discharge was fixed at 9 amps, 
and the axial magnetic field strength was H5 gauss. Their excitation 
mechanism was the same as that employed in the present investigation. 
Essentially the magnetic excitation device employed created sufficient 
variations in the local electron diffusion coefficient to modulate the 
plasma density. (This method of excitation will be explained in detail 
in Section 3.) 

Little and Jones detected the waves by observing variations in the 
light output of the plasma (due to the plasma density variations) by 
means of photomultiplier detection. This method was also used in the 
present investigation and also will be explained in Section 3. Fig. 4 
shows a schematic diagram of the detection apparatus used by Little and 
Jones. 

Their experimental dispersion curves show a dependence of phase 
velocity on frequency and a low frequency cutoff below which propagation 
cannot occur. These results are not directly comparable to those of 
Wong, Motley, and D'Angelo, since the frequency ranges relative to the 


ion cyclotron frequency are entirely different. Fig. 5 shows a plot of 
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FIG. 


с) versus & for 45 gauss external field. The frequency cutoff condition 
can be Seen. 


The results reported by Little and Jones conform to the theory de- 


r 1 


rived in Section 2, although their agreement with their theory is not as 
good as in the present thesis. This is probably due to the unclear 
nature of the boundary conditions in the discharge tube of their experi- 
ment (as they have pointed out). The boundary conditions applicable to 
the present investigation will be discussed in detail in Section 2, and 
it will be shown that they can be unambiguously specified. 

The present investigation examines the dispersion of ion acoustic 
waves in a high density plasma of high degree of ionization (approxi- 
mately 90%) and considers the combined effects of ion drift, damping, 
ion temperature, and a confining magnetic field together with finite beam 
geometry. 

The present thesis has three purposes. First, an improved disper- 
sion relation is derived for a cylindrical plasma including the afore- 
mentioned effects Second, this dispersion relation is tested against 
experiment. Third, from the comparison of theory with experiment, the 
above 5 parameters for the Steady State Plasma Facility at the Naval Post- 
graduate School are computed by least squares computer search to give 
the best fit to the data. The validity of the deduced values of the 
parameters is also demonstrated. 

The derived theoretical solution to the wave equation is a complex 
equation. The real part of the solution predicts the dispersion rela- 
tion, that is wave number as a function of frequency. The imaginary part 
predicts the damping as a function of frequency. 


The solution to the wave equation is derived by assuming the ion and 
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electron pressures are scalars. This allows truncation of the set of 
moments of the Boltzmann equation after the first 2 moments, which are 
the equation of continuity and the momentum transfer equation respective- 
ly. These equations are combined with the plasma equation of state.  As- 
sumption of the azimuthal and longitudinal forms of the wave behavior 
allows the resulting partial differential equation to be reduced to an 
ordinary differential equation involving only the radial behavior. This 
is solved in series form and yields the solution to the dispersion prob- 
lem when the boundary condition is inserted. 

Experimentally the waves are excited by perturbing the local elect- 
ron diffusion coefficient perpendicular to the confining static magnetic 
field. Increasing the local electron cross-field diffusion coefficient 
decreased the local electron density. Since the ions and electrons move 
in phase, this creates a deficiency in the local plasma density of ions 
and electrons. This density perturbation propagates parallel to the 
Static magnetic field. In the present investigation the frequency of 
the applied perturbations was varied from 187 kc to 352 kc. 

The wavelength and amplitude of the ion acoustic waves are measured 
as functions of frequency. The results are reduced by machine computa- 
tion and fitted against theory. As has been previously mentioned the 5 
adjustable parameters in the theory are ion drift velocity, the effective 
ion collision frequency for momentum loss, ion thermal energy, plasma 
beam diameter, and the e-folding distance for the radial density profile. 
The first 3 are characteristics of the plasma facility, while the last 2 
are functions of the applied static magnetic field. Independent deter- 
minations of 4 of these 5 quantities have been made, and these values are 


compared to the results of the least squares adjusted values obtained in 
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the present investigation. The results are presented as functions of the 
static magnetic field, and theories are presented to explain their be- 


havior. 
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Section 2 

BASIC THEORY AND DERIVATION OF THE DISPERSION RELATION 

If Landau damping is to be a factor in wave propagation, the dis- 
persion analysis must begin with the consideration of the particle dis- 
tribution functions, since the mechanism involved is the resonant or 
near-resonant transfer of energy between the wave and those particles 
having thermal speeds approximating the wave phase velocity. If there 
are more particles whose thermal speeds are slightly less the wave phase 
velocity than there are particles whose thermal speeds are slightly 
greater than the wave phase velocity, there will be more particles taking 
energy from the wave than there are giving up energy to it. Then there 
will be a net loss of wave energy, and Landau damping will exist. In 
any system having a Maxwellian distribution there will be some particles 
having thermal speeds approximating the wave phase velocity, but for the 
particular plasma under consideration this is a negligible fraction of 
the total number of particles. 


As will be shown later, the phase velocity of these waves is greater 





than or equal to (нту, Spectroscopic measurements by R. Booth and 
R. L. Kelly of this laboratory have established that ion energies are 
less than 0.5 ev when argon is used in the plasma facility. [10] Pro be 
measurements made by D. M. Gall have established that electron energies 
along the cylindrical beam axis are 22 ev. 11] Thus it may be safely 
assumed that the electron temperature greatly exceeds the ion tempera- 
ture. (This assumption is also consistent with the results of the data 
reduction presented in Section 4.) Thus there will be an insignificant 


number of positive ions capable of causing Landau damping. 


An estimate of the Landau damping may be made by examining a dis- 
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persion relation derived by Fried for the case T>>T. [12] 








m > 
\ : 5 
m 511 1(7 ግ; 
Z AD2 ы : 3 
24( 1 + LE? ας ን s. УЛ” 3/2 
1Y 2 
where 
m is the electron mass 
and 


M is the ion mass. 
For the plasma studied the Debye length is of the order of 0.02 mm, 
whereas typical wavelengths varied from 5 cm to 20 cm. mus £A Av? er 


and may be neglected compared to 1. Making this p— 


4 K 1. ו‎ 
со = ED: (1 = li aedem]: 


Thus it can be seen that an upper limit on the Landau damping term (the 


de 


imaginary term in brackets) is E у 2== 2m δὲ 1073. This is sufficient- 





ly small that it may be ignored compared to the damping arising from ion- 
neutral collisions. This neglect of Landau damping is even more fully 
justified if ion-ion collisions (ion fluid viscosity) are considered. 

The effect would be to introduce disruptive collisions which would des- 
troy the phase relationship between the resonant ions and the wave. This 
critical phase relationship could also be destroyed by ión-neutral colli- 
sions, but not appreciably by ión-electron collisions due to the negli- 
gible electron-to-ion momentum transfer. It should also be noted that 
the thermal speeds of the electrons are much higher than the wave phase 
velocities encountered in the present investigation. Thus the electrons 
see these waves as quasi-static and do not contribute any Landau damp- 
ing. The process of Landau damping therefore need not be considered for 


the purposes of this thesis. The derivation will neglect moments of the 


27 


Boltzmann equation higher than the second. 

For the conditions in this experiment under which ion acoustic waves 
propagate, ions and electrons move in phase due to the strong electro- 
static restoring forces, so the instantaneous ion and electron densities 
are equal. Although the Coulomb interaction is the mechanism by which 
these waves propagate, there is no macroscopic electric field associated 


with the waves. This assumption is justified at any frequency if the 

2 
O 
2M 
(essentially NKT). [15] In this experiment ( 


J), is much greater than the plasma pressure 


2 


magnetic pressure, ( 





) Z 5000 joules/m?, 





while NkT = 10 joules/m?. The following derivation will therefore omit 
consideration of any electric field associated with the waves. 

In the following derivation the quantities W, Во, Т, T are assumed 
to be functions neither of position nor of time. 

W is the equilibrium ion drift velocity 
V is the effective ion collision frequency for momentum loss 
due to ion-neütral collisions 
These assumptions are justified in the cases of By and T. They аге 
probably approximately true in the case of T as well. Since W is defined 
as an equilibrium drift velocity, it must be independent of time, and for 
simplicity it is assumed to be independent of position as well. 

The momentum transfer equation for ions is first linearized and then 
combined with the linearized equation of continuity and the plasma equa- 
tion of state. The usual form is assumed for the azimuthal and longitu- 
dinal forms of the wave solutions, and a series solution of the resulting 
ordinary differential equation of the ion density perturbations is ob- 
tained. Insertion of the radial boundary condition yields the disper- 


sion relation which is then compared to the experimental data. 
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Assuming the electron and ion pressures are scalar quantities, the 
momentum transfer equation for the ions under the influence of the per- 


turbation is: 





x Тара EG.» NJ | азоб +] (5) 


+ (N+ п) (У + WW “IN + п) (У + אכ‎ 


where 


is the equilibrium electron pressure 


Ay 


is the equilibrium ion pressure 


is the electron pressure perturbation 


A ю 


is the ion pressure perturbation 


<| 


is the ion velocity induced by the wave 
n is the ion density perturbation 
and the positive 2 direction is taken in the direction of Bo- 
The momentum transfer equation for the ions in the unperturbed 


plasma column is: 
zl ge JD -| 2— + w q ][ 6] 5 м2. оме (5) 


Expanding equations (H) and (5), subtracting (5) from (H), and 
assuming the density perturbation is small compared to the unperturbed 
density allows the resulting equation to be partially linearized to the 


first order in small quantities. 
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Since solutions of the form гі 9 +2 гар are sought, the follow- 


ing operators may be replaced by the listed variables: 


ዊ а -ic) (c9 assumed real) 


2 


id (4 complex longitudinal wavenumber)‏ چس 


2 
1 complex azimuthal wavenumber 
SS ger ο ad 
ЕЕЕ ^ 2 
Also since W is assumed to be in the z direction, W - WZ and 
vou. 
w-V SS 
The radial form of the density perturbations is left unspecified for the 
Present. 
Linearizing the expression V. (М + п) (У + W to the first order 
in small quantities gives: 
== — -- e es אכ‎ 1 . Ë: 4 
у “У (N+ ny (V + W = W ርኛ መ ነ a nVg * inv», 


where 


^ 
Ver + Vg8 + V,2. 


Thus (6) becomes: 


ም - ON 
-Vp +w) = -iONV - iodnW + WV, ST 2 (7) 
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From (7) are immediately obtained: 
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Forming the vector difference of equations (7a) and (7b) and taking 


the divergence of the resulting quantity gives: 


Әр 2929, 1,3. d - 
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Neglecting volume ionization and recombination the equation of 


continuity for the ions is: 
й 2-ə---2 ? 
. (N + n)(V + W) = - sr NN +0) = ісп (9) 


Taking the divergence of equation (7) gives: 


=s +) = 7 NV + Y - ag dani: 2 Ur i 


(10) 
+ idy[D nT eT. a] [v.e vens] 
- | 86 - ነ" | 
Linearizing equation (9) gives: 
V. (N+ n (V + WÀ = Ç. NV + Ὁ. τν (11) 


Using equation (11) to evaluate V. (NV) + ў. (NW 


е е А 
and using equation (8) to evaluate V . (NVgr - NV.8), 


equation (10) becomes: 
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Under the assumption of the spatial invariance of Y and T equation 


(12) becomes: 


n 1 Dn Ма (14) - if 226 2 
(13) 


+ (с^ + с? 


шт iL- (і 448 - 2) + 4 
א‎ Or Ga- ilu- α΄ 
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where c^ a — (definition) 


and C2 = - (definition). 


The differential equation (13) would be directly reducible to 


Bessel's equation of order M if were equal to zero (that is if the 





ON 
or 
unperturbed column had a uniform radial density distribution). This, 
however, is not the case in the present investigation. 

The mean free path for the ion-neutral collision process is greater 
than l m, while the ion gyro radius is less than 1073 m. Up to several 
ion gyro radii away from the beam axis the analysis of Pfirsch and 
Biermann predicts that for a Maxwellian velocity distribution and a 


Gaussian distribution of the guiding centers of the charged particles 
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there will be a Gaussian density profile. [14] Previous probe analysis 
by Gall and Oleson has fitted a Gaussian radial density profile to N, 


22 
at least to a good approximation. [15] Thus א‎ = Noe `Š Б » so that 


> SE = As and equation (13) becomes: 
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equation (14) becomes: 
22 > 352 р 
52 ІП... «ааа Шы. ар d ыы + n(s? - 2) = 0 (15) 
ο we G3-QW * ip) ай 


Assuming a solution of the form 
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and inserting this into equation (15) leads to the recursion relation 
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Thus the solution of equation (15) is: 
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Equation (16) reduces to the Bessel function Jn for the case g = 0 
(uniform density in the unperturbed column). For this case the disper- 
sion relation would assume the simple form \X = 2.405 where the boundary 
condition has been imposed that r = £ is the radial position at which the 
perturbation vanishes. The validity of this type of boundary condition 
was confirmed by the probe measurement discussed in detail in Section 3. 
Independent confirmation that this type of boundary condition is appro- 
priate can be obtained from an examination of the radial electron temper- 
ature profile as obtained by Gall and Oleson. 15) This curve is repro- 
duced as Fig. 18 in Section 5. It is evident from an examination of 
Fig. 18 that the electron temperature drops sharply (or at least becomes 
small compared to the value at r = 0) at the position r = 7,8 mm. Ап 
examination of equation (14) under the condition that both c and C be- 
come small reveals that the magnitude of the density perturbation, n, 
must also approach zero. The use of this type of boundary condition is 
particularly convenient, since no assumption need be made concerning 


the existence of either a pressure node or a velocity node at the plasma 
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boundary. This problem has plagued investigators working with low ion- 
ization gaseous discharges. Б 


For the case ₪ < 0 the dispersion relation for /(- O is: 


со 1 | 
>. ቦ (1 - “ን (አሮ)? 


= = 0 (17) 
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where now 
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Equation (17), or more accurately the recursion relation between 


and 


"א 


the coefficients, is compared to the experimental data by a least 


Squares direct search computer program on a CDC 1604. 
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Section 3 
EXPERIMENTAL ARRANGEMENT 

A schematic illustration of the Steady State Plasma Facility at the 
Naval Postgraduate School appears in Fig. 6. The plasma is produced by 
a hollow cathode discharge, which has been investigated by Lidsky et. al. 
[16] The details of its construction appear in Fig. 7. The facility is 
capable of operation on helium, neon, and argon, although only argon was 
employed in the present experiment. The neutral argon gas is supplied 
through a micrometer valve assembly to the hollow cathode. An r-f volt- 
age is supplied between the cathode and the first anode in order to ig- 
nite the arc. Once the arc has been struck, the discharge is self- 
sustaining, and the r-f voltage is turned off. 

The plasma so produced drifts through the hole in the first anode 
and into the experimental tube. This section is a 3 meter pyrex tube of 
inside diameter 10 cm. ፲ቲ is axially aligned with the confining magnetic 
channel to within 1 mm throughout the length of the experimental region. 

The optimum operating conditions for the facility are an arc cur- 
rent of 100 amps and a voltage of 90 volts between the cathode and first 
anode. This condition is most conveniently achieved by regulation of 
the argon g£as pressure through the micrometer valve assembly. 

Differential pumping of neutral gas atoms by the diffusion pumps 
maintains a neutral gas pressure of approximately 3 x 107? mm Hg through- 
out the experimental region. 

The strength of the confining magnetic field is continuously variable 
from 500 gauss to 10,000 gauss. Some experimental data were taken at a 
field strength of 600 gauss, but the agreement with theory was quite poor. 
It was suspected that some drift occurred in the current supplying the 


field coils, so that the confining magnetic field was not constant with 
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time (a necessary condition for taking a complete set of dispersion data). 
These data are not reported, and the minimum field strength employed was 
800 gauss. 

Upon reaching the end of the experimental region, the plasma ions 
and electrons strike the floating anode (so-called because its potential 
is allowed to float) and are recombined. They are then pumped away as 
neutrals by the vacuum pumps in the burial chamber. A photograph of this 
section appears in Fig. 8. A side view of the entire facility appears 
in Fig. 9, and a view of the cathode-first anode region appears in Fig. 
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In order to excite the wave, a short r-f coil was tightly wrapped 
around the pyrex tube as shown in Fig. 6. It was designed for optimum 
coupling to the plasma at the center of the band of exciting frequencies 
employed. This coupling effectiveness was determined by the length of 
the coil, which was approximately that of the wavelengths excited.  Ex- 
perimentally it was found that the r-f field strength dropped appreciably 
in the last % of the coil at each end, so that these outer regions were 
ineffective in producing a perturbing field. The coil itself was approx- 
imately 15 cm long, leaving approximately 7 cm for the effective field- 
Producing region. This was in the center of the band of wavelengths 
excited. 

The coil was the inductive arm of a parallel "almost resonant" tank 
circuit, which was tuned (with approximately 1500 мжЕ excess capacity) 
to resonance at the applied transmitter frequency. This slight departure 
from the resonance condition was necessary so that the transmitter would 
See a load impedance similar to that of the radio transmission antenna 
which it was designed to feed. 

The transmitter used to power the tank circuit was a U.S. Navy 
type TAB-7 manufactured by the Westinghouse Corp. It is capable of 
producing a peak r-f power level of 2 kw over the frequency range of 
100 ke to 500 kc. This full power condition was not used due to the 
inability of the r-f exciting coil to dissipate the heat generated at 
this condition. Typical operating conditions were a transmitter output 
current of 3-4 amps. The Q value of the tank circuit was estimated to 
be about 10 at the frequencies used. Thus the actual exciting current 
in the coil was approximately 30-40 amps. Direct measurement of the 


r-f magnetic fields produced yielded values of from 125 to 195 gauss, 
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depending on transmitter power. As previously mentioned, this perturbing 
magnetic field produced variations in the local electron diffusion co- 
efficient, creating plasma density variations. 

Assuming the magnetic perturbation, B, is of the form B - B,sin WT 
and that J). (which was true in all cases), only the influence of the 
perturbing field on the electrons need be considered, since the massive 
ions are unable to respond directly to the rapidly varying magnetic field, 
while the electrons can instantaneously follow its variations. The dif- 
fusion rate of electrons perpendicular to the magnetic field is: 


Dp = 


1 + [ SCH + ye ir 


where Dp is the local electron diffusion coefficient with both B and 
Bg applied 
D is the local electron diffusion coefficient with zero applied 
magnetic field 
e is the electron charge 
m is the electron mass 


N is the electron-neutral collision frequency 
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a is an average effective diffusion mean free path, there results 


αμ or N = RI . Thus N& 
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sin dF, and the plasma den- 





sity is modulated at the applied frequency, the strength of the perturba- 
tion being proportional to that of the r-f magnetic field. These per- 
turbations were studied over the frequency range of 187 kc to 352 kc and 
for static magnetic fields from 800 to 2200 gauss. 

Once the waves have been generated, they propagate parallel to the 
Static magnetic field in both directions. Only propagation anti-parallel 
to the static magnetic field was investigated in detail. (One set of 
parallel propagation measurements was made; this allowed an independent 
determination of the ion drift velocity, W.) 

Previous spectroscopic measurements have established that virtually 
all the light emitted by the plasma is from excited neutral atoms. [10] 
Therefore the emitted light intensity at any point along the plasma col- 
unn is directly dependent on the local electron density. (An ion- 
electron light production mechanism would produce a light intensity which 
was proportíónal to the square of the local charged particle density.) 
Thus a photomultiplier tube can be used to detect the presence of the 
waves. The detection circuit schematic is given in Fig. 11. The photo- 
multiplier tube was housed in a light-tight electrostatic and magnetic 
shield manufactured by the Quanta Corp. A 1P21 photomultiplier tube was 
mounted inside this shield and dynode voltages were supplied by batteries. 
The entire assembly was located in a region of minimum magnetic field 
(less than 25 gauss). The optical signal was fed to the phototube 
through a 10 ft flexible unaligned optical fiber bundle manufactured by 
Optics Technology, Inc. This fiber bundle caused a 67% loss in light 


intensity, but the signal strength was large enough that this caused no 


45 


SSS 
XINOYLY31 


O 


ፎ0ዬ qp O ¿+ 
INS VOS * dH 


ዕ11 
AS IHLISN 00 +ל0‎ 








2 - 1 
310NN8 Y38l3 1VOIl1dO 3181X313 


— 300NV 
ברסצן ואפ‎ 


3003142 >--- 





318ΥΛΟΝ 


JILVN3HOS כוטכחו1‎ 2 
ІІ 914 





problem. 

A Keithley Model 110 Wideband Amplifier with an overall gain of 
unity was used as an impedance matching device between the phototube 
and the Hewlett-Packard Model 450A Wideband Amplifier (which was set to 
a 20 db gain). In order to improve discrimination against noise the 
signal was then fed into a Spencer-Kennedy Electronic Filter Model 302, 
which served as a band pass filter centered at the transmitter frequency. 

The next element in the detector circuit was a Hewlett-Packard Model 
400H Vacuum Tube Voltmeter which was used to make amplitude (attenuation) 
measurements. Finally a Tektronix Dual Beam Oscilloscope Type 555 was 
used to display the detected signal and compare it with a monátor signal 
taken directly from the TAB-7 transmitter. Both signals were subject to 
various phase delays (due to coupling into and out of the plasma, ele- 
ments in the detector circuit, etc.), but comparison of the phases of 
the two signals appearing on the oscilloscope gave the relative phase 
delay of the arriving signal. Moving the optical fiber bundle along the 
plasma column until the phase of the arriving signal was the same as be- 
fore allowed a direct determination of the wavelength of the propagating 
signal. 

Signal strength after amplification and filtering was of the order 
of 100 mv (depending on dynode voltages and transmitter power).  Signal- 
to-noise ratios varied widely depending primarily on the electronic 
filter settings and on the transmitter power. Individual filter settings 
at each frequency employed were necessary in order to discriminate 
against the noise background in the plasma. The frequencies selected 
for study were chosen in the quieter regions of the plasma noise spec- 
trum. There was a marked variation in the location of these regions 


depending on the static magnetic field strength. A trial and error 
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approach proved to be the most productive method for finding their loca- 
tion. 

A typical oscilloscope photograph is shown in Fig. 12. The upper 
trace shows the detected signal (after amplification and filtering). 
Part of its amplitude is due to an unavoidable inductive pickup condi- 
tion. This part was, of course, subtracted from the raw data before 
analysis. Some plasma noise can be clearly seen within the signal en- 
velope. The signal-to-noise ratio here is slightly greater than 1. The 
lower trace is a comparison signal taken from the TAB-7 transmitter. 
This photograph was taken for a 250 kc signal. Deflection sensitivity 
of both traces was 100 mw/cm, and sweep speed was 5 Msec/cm. 

Wavelength and amplitude measurements were taken as functions of 
frequency for various magnetic field strengths. The results provide 
the raw data from which the dispersion and attenuation curves presented 
in the next section were derived. 

To insure that modes other than the y= 0 were not excited, probe 
measurements were made in an attempt to detect the presence of the 


ΛΑΞ 1, 2, 3 modes. (The form of these modes is sketched below.) 


ΑΞ! M=2 ሥ=3 
At maximum detection sensitivity no such waves were detected. Since the 
excitation mechanism had azimuthal symmetry, modes other than = 0 


should not be present. Actual confirmation that the detected wave was 


actually the u = O mode was made by probe measurement. The perturbation 
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Figure 12. Typical Oscilloscope Trace 
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was observed to have a maximum at r = 0 (on the beam axis) and was ob- 
served to decrease as the probe was moved out from the axis. This be- 
havior would only be possible for the /^7 O mode. All other modes would 
exhibit a null at r = 0. This probe measurement also forms the basis for 
the comparison theory for the cutoff distance (that is Í). These mea- 
surements showed no detectable perturbations beyond a distance of 15 mm 
from the beam axis. 

No standing waves were observed to exist. This was probably due 
to the effect of the burial chamber surrounding the floating anode. In 
the burial chamber ions and electrons are recombined and the resulting 
neutral gas pumped away. Thus the plasma charge density is low, and 
the neutral gas density is high. Both of these conditions would tend 
to preclude the establishment of a standing wave pattern. Reflection 
of the waves from the floating anode could produce a standing wave pat- 


tern, but this effect has never been observed. 17 
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Section Y 
EXPERIMENTAL RESULIS 
Given a fixed static magnetic field strength, a set of data consists 
of @ and ሥ as functions of eJ. The raw data are reduced in the follow- 
ing way. Assuming © is real, let y^ = OX + 1 6 ( & and 6 both real). 
First the measured wavelengths are reduced to dimensionless variables 
«c 


by defining y z ER Knowing the wavelength allows direct calculation 





of the damping decrement , (6 . If two intensity measurements are made 








I x E E 
one wavelength apart, l = e ፆ|።1 28. |21 - 2 - T, Thus 
Io ο ol 
I _ ጋዝሪ 
= =e: Y , from which f can be calculated directly. The damping 
O 


decrement is then reduced to a dimensionless variable by defining 


а = 4e . The frequency is reduced to a dimensionless variable by de- 


€) 


fining x = >. For a given magnetic field strength the entire set of 


data (that is y and z vs. x) are read into a direct search computation 
program. The 5 adjustable parameters W, V, c c and Š are reduced to 


dimensionless variables by defining 


c Ë. Š c 
Eu coge 
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The computation program símultaneously adjusts these parameters to 
find the dispersion function yielding the least square error when com- 
pared to the data. In addition to performing this calculation, graph- 
ical displays are also presented for y vs. x and z vs. x. The plotted 
Squares represent the data points, and the continuous curve represents 
the dispersion relation using the 5 parameters just computed. Results 
are presented from 800 to 2200 gauss in the following curves, the para- 


meter being Bo (web/m2?). 
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The following 4 curves present the bahavior of Y (the dimension- 
less drift velocity) and @ (the dimensionless ion thermal speed) plus 
the non-normalized parameters W and C as functions of the magnetic field 
Strength. Independent measurement of W was made at 1000 gauss by study- 
ing parallel propagation, and the results of the data reduction for that 
set of dispersion data gave a value within 12% of that obtained from the 
anti-parallel analysis. This will be discussed more thoroughly in Sec- 
tion 5. However, no theories are presented to explain the behavior of 
either parameter. A program was prepared to fit a least squares straight 
line to the deduced parameters as a first approximation. These lines 
are the continuous curves, and the squares represent the fitted para- 
meters for each magnetic field strength. The standard deviations of 
both slope and intercept for these plots appear in Table 1. The units 


of each of the non-normalized quantities appear on the individual curves. 


Table 1 
C ion 
Quant ity 7 = - gQ” = === W(m/s) energy (ev) 
std. dev. of slope .00168 .0741 12:2 „Ш20 
Std. dev. of intercept .000264 . 0116 1.9 .0188 
slope . 00109 .415 7.92 27 
intercept -.0124 .0614 -90.0 ‚00635 


59 


914 


212 


ГА 


BENI 
ΜΡ ΠΡΙΝ 
ιτ 
таыпттны 
ИІ 
ከከ | 


8 


Вас 


E 906 


O 
@ 
с 
po 
со 
O 
፦ዶ 
σι 
ኢን 
O 
ኤን 
31 


N SCRLE = 5.66፻ ብ2 ህከ113“ጉርዚ 
Y-9CRLE * 2.00E-02 LINITS/INCH. 


STEMA ÛU STATIC MaGNF TTC. FTFI T) CAIER^ZSQ ME T) 


60 


K-SCRLE - B.00E-02 UNITS/INCH. 
Y-SCALE = 2.00E--01 UNITS/INCH. 


NF DRIETOUEPBOCITY (uz Y 
STATIC MAGNETIC FIELD (WEB-SQ MET) 


61 


| y 
ስ 
noo 865 010 916 020 R25 839 825 


K-SCRLE = b.00E-02 LNITS^INCH. 
Y-SCALE - b.00E-02 UNITS~INCH. 


ION THERMAL ENERGY (EU) U 
SE elle Zee (ЕВ 54 МЕТ) 


X-SCALE = 5.88E -82 Pos d. 
Y-SCALE = 2.00E -02 


VU STATIC MAGNETIC FIELD (WEB-3Q MET)  *‏ א 


The next set of 6 curves represents the behavior of 5 (the dimen- 
sionless effective ion collision frequency for momentum loss), & (the 
dimensionless radial cutoff distance for the perturbations), and € (the 
dimensionless e-folding parameter for the unperturbed radial plasma den- 
sity distribution) plus the corresponding non-normalized parameters Y, 
Ў, апа б. Again the squares represent the data points deduced from the 
data reduction, and the continuous curves represent the theories for the 


behavior of these parameters. These theories are presented in detail in 


Section 5. 
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The following set of curves depects the behavior of the dimension- 


less phase velocity, 22 ‚ апа the dimensionless group velocity, 


Els » for each value of static magnetic field strength. Again the 
parameter is the value of B, (web/m2). The squares represent the normal- 
ized data points and the continuous curves represent the values of S 
and £g as taken from the dispersion relation using the values of the 
5 parameters appropriate to the particular value of Bg. 

The modulation feature of the TAB-7 transmitter was inoperative, 
So no group velocity measurements were made, and there is no data with 


which to test the theory. For this reason no data points appear on the 


group velocity curves. 
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Ihe last ^ curves represent the general behavior of the dispersion 


function throughout the first quadrant of Á , espace. Plots are pre- 


voc = ይ. x3 — + vse x, and — 5 





sented for y = ха 
vs. x. The 1000 gauss magnetic field setting was typical, and these 
plots represent the dispersion function with the 5 parameters having 
values appropriate to that setting. Since the 1000 gauss data points 
have already been presented in the previous curves, the data are omitted 
from these 4 curves. 

The first 3 curves are for the complete quadrant, but the group 
velocity is depicted only for the regions in which the damping is small, 
that is where ion acoustic waves might be useful for communication pur- 
poses. These are also the regions where 227 O. In the regions of 
anomalous dispersion (Se <0) the concept of group velocity becomes 


vague, since it no longer necessarily represents the velocity of propa- 


gation of a signal. [9] 
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Section 5 

INTERPRETATION OF EXPERIMENTAL RESULTS 

Before proceeding to a more detailed analysis of the experimental 
results, the basic behavior of the dispersion curves and the attenuation 
curves will be discussed. Consider first Fig. 13; here a non-dimension- 
less plot of X vs. GD is presented. For simplicity, and in order to 
isolate the effects of the various parameters on the basic dispersion 
curve, consider first the case of W = 0 (zero ion drift velocity); 
2: 0 (по ion-neutral collisions and A = С + ip is real); C = O (cold 
ions); and б = O (uniform radial density distribution ig the unperturbed 
plasma column); the dispersion relation, equation (17) reduces to: 


02 2 су? - Xo? 
O ммм 


where ጻዕ ₪5 the firsi$ root of Jo x == 
This equation is plotted in Fig. 13 for the 4 distinct cases: 
41720, 4 = © (zero magnetic field, infinite medium) 
270, R=0 (non-zero magnetic field, infinite medium) 
С. = 0, ያ # 20 (zero magnetic field, finite medium) 
LAO, LAO (non-zero magnetic field, finite medium) 

The cases presented in Fig. 13 are elementary and are presented 
as an introductory clarification for the more detailed analysis which 
follows. 

From Fig. 13 it is evident that the presence or absence of a mag- 
netic field has no effect on the dispersion in an infinite medium. The 
reason is that the perturbational velocities Vg and V, are always zero 


for this case (see equations 4a and 4b). Only V, is non-zero, and since 


it is not a function of sL (see equation lc), there is no dispersion. 
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The term cutoff will be used to describe the propagation condition 
under which the phase velocity of the wave is infinite (or when the real 
part of the propagation constant, A , is zero). The term resonance will 
be used to describe the condition under which the phase velocity is zero 
(A is infinite). 

For a finite medium the high frequencies (and hence small wave- 
lengths) can easily propagate down the plasma column. As the wavelength 
increases there is a point at which the plasma can no longer transmit 
the signal. This appears as a cutoff at frequency 621 - መ፦- , апа 
means that at frequencies less than €, «? «0 meaning propagation is 
impossible. This is a consequence of the location of the boundary sur- 
face (the cylinder of radius f Dis 

Finally for the case of a finite medium with a non-zero magnetic 
field, there is a resonance at the ion cyclotron frequency. At this 
frequency all the ion momemtum appears as cyclotron motion. Again 
e x At fre- 

2 


quencies ALALA, the plasma beam cannot transmit the signals, and prop- 


there is a cutoff, but this time at D AT (92 + 


agation is impossible since ©? < 0. 

Tn the following discussion it will be assumed that the dispersion 
relation is exactly the Bessel function of order zero. This is equiva- 
lent to the assumption that the unperturbed plasma column has a uniform 
radial density profile. This will considerably simplify the discussion, 
since the infinite series for the dispersion relation need not be consid- 
ered. For this case the relation እሪ = 2.405 holds. In actual practice 
the perturbations differ from pure Bessel functions by no more than about 
10%, so the approximation is justified for discussion purposes. 


From this point on the assumption will be made that the plasma is 
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finite in the radial direction (XZ o0) and that there is a non-zero 
magnetic field (41 X 0). Also all variables will be normalized (in- 
cluding the 5 parameters) to the dimensionless values. These definitions 


are repeated here for clarification. 





= с _ ON 
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The effect o (lion drift) is easily understood, since it intro- 
duces a Doppler shift due to the moving medium. In this experiment 7 
and M were anti-parallel, so the phase velocities are slightly reduced 
by the drift. This is illustrated in Fig. 14. 

A set of measurements was made with 7 and 4 parallel at a field 
strength of 1000 gauss. Direct measurement of the difference in phase 
velocities for the parallel and anti-parallel cases (not the use of com- 
puter data reduction) gave an entirely independent determination of the 
value of the ion drift velocity. This value was within 12% of that ob- 
tained by the computer data reduction process, and this agreement lends 
strong confirmation to the values indirectly deduced from the measure- 
ments. 

The effect of ion temperature is also easily understood. Since ion 


acoustic waves are density perturbations, they propagate at a velocity 


(Pressure term ₪ 


: š Passing from the case of cold ions 
inertia term 


proportional to 
to warm ions will increase the pressure term by a small amount equal to 
the ion pressure. Thus the phase velocity will be higher for a given 


frequency; that is 4 will be reduced. This effect is illustrated in 


Fig. 15. The solid curve represents the case of cold ions, and the 
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dashed curve represents the case of warm ions. 

The effect of { is more involved. If ion-neutral collisions are 
considered, there will be a decrease in the wave phase velocity to a value 
smaller than in the collisionless case. [7ב]‎ In the present investigation, 
however, this effect is of the second order. The more important effect 
of collisions, however, is that a damping mechanism has been introduced. 
This means that there is no longer an undamped resonance at x ። 1. Rather 
the resonance is diffuse, and the maximum wavenumber is limited due to 
the damping. Waves can also exist in the previously forbidden region, 

1 > x < Ses ; although here they are heavily damped. This effect is 


illustrated in Fig. 16. Although attempts to excite and detect waves in 





the region 14 x 494 were made, the heavy damping plus the long wave- 


lengths in this region precluded their detection. 
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Топ acoustic waves are produced by the Coulomb interaction (a long 
range interaction). Propagation is possible only when the ion-neutral 
collision frequency is small compared to the wave frequency. The few 
collisions present disrupt the long-range ordering, and damping increases 
with increasing collision frequency. Typical ion-neutral collision fre- 
quencies for this investigation were less than 500 бесті, so the conti- 
tion of heavy damping would be encountered only at frequencies much 


less than those studied (plus the heavily damped waves in the region 





כ6 

| < x = ) 9 

Due to the apparent invariance of 7 with magnatic field, it is sus- 
pected that the mechanism producing the drift is an electric field in the 
interior of the plasma, rather than a fluid pressure differential between 
anodes. If the ion drift were due to the latter mechanism, its value 
would depend on the magnetic field strength, since the neutral gas in- 
jection pressure varied with the field strength. Unfortunately previous 
electric field measurements were not available, so a confirmation of the 


mechanism which produces the drift cannot be made. 


“2, 22 
For ion-neutral collisions, A. = 1 суш... 
ion ING כ‎ 
8 
where א‎ is the ion mean free path 
Ng is the neutral gas density 
@ is the collision cross section for the process 
Vo is the mean ion thermal speed 


Knowing σ-. and Ns (which is manually adjusted for each magnetic 
field strength to give optimum machine operation), a straightforward 


calculation a D - O7 Pg x 246 x 10° where Ες is the neutral gas 


ν' 
ς 


е 


pressure and Ў” = 


ЭМ! 


5 ν) 


Bo 





Thus = 41.7 x 10 


In the above calculation it is assumed that there is a 100% excess ion 
momentum loss per collision. This assumption is not strictly true, but 
it is certainly reasonable, since a large fraction of the directed mo- 
mentum is lost in any collision. Tt is also assumed that the ion per- 
turbation velocities are small compared to the mean ion thermal speed. 
Thus the collision frequencies will not depend on the strength of the 
perturbation. 

The above equation appears as the theory on the ion temperature 
curves in Section 4. The value of Q used was taken from the data re- 
duction program appropriate to the magnetic field strength. 

Since TK 1 (the ions are much colder than the electrons), the dis- 
persion equation is quite insensitive to changes inU . This can be 
seen by examining the data points of the U^ curves in Section 4. It is 
considered by the author that the deduced values of V^ (and hence the 
ion enérgies) should not be interpreted too literally. Іп the V^ vs. BS 
curve the deduced data points are rather widely scattered. This inac- 
curacy is reflected not only іп the V" vs. B curves, but also in the 
₪ VS. Bg curves, since J x0. But it can be seen, at least qualita- 
tively, the observed damping results primarily from ion-neutral colli- 
sions. 

Variations in Bj affect the propagation characteristics (that is 


22 
(Э, = (0? + x ) 2, since £L is of the same order of magnitude as 


ር 
na An analysis by Wong [10] concludes that the transverse boundary 
conditions are important only if the ion cyclotron radius, t is com- 
parable to the plasma beam radius. At 1000 gauss f^i = 3.2 mm 2888 = 15 


ша, so this criterion is approximately met. 


ЭА 


Neglecting % (uniform unperturbed column) Wong has shown that 
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where Z' is the derivative of the plasma dispersion function 
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і : አጅ тэ % 055 | 
In the present investigation, 4 > 26 2 
-. γα” х02012 


Thus Z' EP = 6 א‎ 0 | and the effect of the radial boundary 
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is important only if n а . Since this is indeed the case, 
X o i? 
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This is a lower limit for the ion temperature, since any tempera- 


it must be true that T 2 


ture appreciably lower would predict that 62 כ‎ is not a function of Bo- 
This estimate is considerably higher than the results of the data re- 
duction (.05 ev to .17 ev), and the spectroscopic measurements of Booth 
and Kelly. 10] For that reason it must be stated that the ion energy 
cannot be determined accurately by this experimental method. It is 
probable, however, that the ion energies lie somewhere between .05 ev 
and 2 ev. The ion temperature has been found to increase with distance 
from the plasma source in a helium plasma, varying in value from about 
l ev at the first anode to about 8 ev several meters downstream. ₪ 

The increase of ion energy with increasing magnetic field is be- 
lieved to be due to the increased electron-ion collision frequency in 
the denser beams resulting from higher field strengths. 

It should be noted that volume recombination and ionization would 


increase the damping by producing hot neutrals and cold ions respectively. 
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Also, since most ion momentum is directed along the К ахіѕ, е1аѕііс 
Scattering of an ion through anything other than a small angle would 
appear as a damping effect. Charge exchange would substitute a hot 
neutral and a cold ion for the original hot ion and cold neutral. Al- 
though the mean free paths for charge exchange and ion-neütral collisions 
are of the same order of magnitude, the actual collision rate for charge 
exchange would be somewhat less than that for elastic collisions, since 
not every ion-neutral collision would result in charge exchange. But 
charge exchange would certainly appear as another damping mechanism. 

All these processes probably contribute to the inaccuracies of the above 
theory. 

It should be noted that the observed damping is much less than that 
observed by Little and Jones in a gaseous discharge where the percentage 
ionization is low and there are many more neutrals present. B Also in 
such a discharge collisions between the ions and the container walls 
would cause additional damping. Such an effect is not a factor in the 
present investigation. 

The primary effect of € is to alter the form of the radial per- 
turbation from J, to a more rapidly decreasing function of r. The effect 
of € on y and z is a complicated second order process; it will not be 
discussed further here. 

Neglecting the second-order effects (7, V , €), the equations 


Sb, and 4e become (for M = 0): 
2p 
=<ራ — 
Gay co 
2р (14 -p) 
A 
(1ὁ- (2 +2 


MNVg = 
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ММУ, = ikp 
= (162 - 9) 


In this form the equations clearly illustrate the damping effect of V. 


If Y is now set to zero, these equations can be solved to give: 
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Since = « Go (except near cutoff), it is clear that most of the 
perturbational momentum appears in the motion of the ions in the В аігес- 
tion. If is for this reason that the waves are called acoustic waves. 

Theory for the € , ሪ curves is derived from the probe analysis 
mentioned in Section 3. Those measurements (which confirmed the pre- 
sence of only the /А - 0 mode) showed that no density perturbation was 
detectable beyond 15 mm from the plasma beam axis. This distance was 
also found to be independent of the magnetic field strength up to 2200 
gauss. This was a remarkably strong confirmation of the values of 
deduced from the data reduction. It also demonstrates the validity of 
the basic dispersion relation, since this relation is far more sensitive 
to small changes in Á than for any other parameter. 

Fig. 17 is reproduced from the work of Gall and Oleson [15]; it can 
be seen that a value of 4 = 10 mm would be a more reasonable estimate, 
at least from the point of view of electron temperature alone. However 
the value of P« - 15 mm taken from the probe analysis gave a much better 


fit to the data and was used as the theoretical value, since it was 
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independently confirmed. 

Theory for the € , 5 curves was taken from the work of Gall and 
Oleson. (15 Although their work did not directly yield values of > as 
a function of the magnetic field strength, it did indicate that this 
quantity is not a sensitive function of the magnetic field strength in 
the range used in the present investigation (800 to 2200 gauss). Fig. 
18 is reproduced to illustrate this point. From Fig. 18 § is seen to 


be equal to 100 m! at pe 600 gauss. Assuming $ is not a function 


230 
of Bo, € Li  "' 
0 


Table 2 summarizes the results of the present investigation for the 


low, medium, and high field strengths used. 


Table 2 

quantity 800 gauss 1600 gauss 2200 gauss 
wavelength (m) -0458 to 1-19 .0314 to .341 .0248 to .145 
i/e attenuation dist. 20 to 500 133 to 1110 455 to 2500 
(m) 
phase velocity/c m 25 to 26 1.25 to 9.0 1.2 to 4.4 
group velocity/c 209 to p62 sl to TAS EI tomo 
ion drift velocity 90.8 89.4 1 
(ms 1) 
effective ion collision 34 2 510 431 
frequency for momentum 
loss (sl) 
ion energy (ev) ‚05 1.3 16 
radial cutoff 203152 . 0147 . 0152 
distance (m) 
radial 1/e profile . 0100 ‚0105 „0105 


parameter (m=!) 
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Since the values of y and z vs. x cover such wide ranges, a tabula- 
tion of these quantities vs. x is not included at this point.  Represen- 


tative values can be read from the curves presented in Section A. 


Table 3 
Parameter Measured Value Deduced Value 
| ל‎ | 95.2 m/s 85.0 m/s 
Y 421 1 480 s^! 
T none .092 ev 
15 mm 15.5 mm 
Š 100 m^! 102 ₪1 
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Section 6 
CONCLUSIONS 

In conclusion it may be stated that for the propagation of ion 
acoustic waves where Landau damping is not present the new dispersion 
relation derived in Section 2 is capable of predicting the propagation 
characteristics of these waves quite accurately. The greatest error 
present in the dispersion and attenuation curves (that is the greatest 
departure between data and theory) is 6%. 1% was, however, much more 
difficult to arrive at satisfactory theories for the 5 adjustable para- 
meters. The close agreement in the dispersion and attenuation curves, 
however, is believed to result from the unambiguous boundary condition 
Specified and from the use of photomultiplier detection (rather than 
probe detection with its inherent errors). 

Where damping processes other than ion-neutral collisions occur, 
the effective ion collision rate for momentum loss may be redefined to 
include all such processes (Such as the collision rate of ions with the 
walls of the container). For this case Y would be the sum of the in- 
dividual collision rates for all processes. Since negligible momentum 
is carried by the electrons, electron momentum loss mechanisms may be 
ignored without appreciable error. 

The externally accomplished excitation and dispersion measurements 
described have been a valuable diagnostic tool for the investigation of 
the plasma characteristics. It was by this method that the values of 
w, 2 7 2 апа 5 were deduced from the data. This method could be 
even more valuable where internal measurements (such as probe diagnos- 
tics) are impractical. 


It is recommended that future investigations be made to test the 
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theory at very low frequencies (1 ke to 100 kc). Due to the frequency 
limitations of the TAB-7 transmitter employed, measurements in this 

lower frequency region were not possible. Future investigations in this 
region would provide an additional test of the theory and additional con- 
firmation of the deduced values of the 5 parameters considered. 

It is also recommended that the noise spectrum of the plasma be 
analyzed. This would próvide additional information on the naturally 
occurring interior phenomena. 

Investigation of ion acoustic wave propagation in other rare gases 
(helium and neon for example) would also provide additional checks of 


the theory. 
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